arXiv:1503.06486v2 [math.PR] 9 Mar 2016 


A fractional counting process and its 
connection with the Poisson process 

Antonio Di Crescenzo, Barbara Martinucci, Alessandra Meoli 

Dipartimento di Matematica, Universita degli Studi di Salerno 
Via Giovanni Paolo II, n. 132 - 1-84084 Fisciano (SA), Italy 
Email; {adicrescenzo, bmartinucci, ameoli}@unisa.it 

To appear in ALEA Latin American Journal of Probability and 

Mathematical Statistics 


Abstract 

We consider a fractional counting process with jumps of amplitude with 

k G N, whose probabilities satisfy a suitable system of fractional difference-differential 
equations. We obtain the moment generating function and the probability law of the 
resulting process in terms of generalized Mittag-Leffler functions. We also discuss two 
equivalent representations both in terms of a compound fractional Poisson process and of 
a subordinator governed by a suitable fractional Cauchy problem. The first occurrence 
time of a jump of fixed amplitude is proved to have the same distribution as the waiting 
time of the first event of a classical fractional Poisson process, this extending a well-known 
property of the Poisson process. When k = 2 we also express the distribution of the first 
passage time of the fractional counting process in an integral form. Finally, we show that 
the ratios given by the powers of the fractional Poisson process and of the counting process 
over their means tend to 1 in probability. 

Key words: Fractional difference-differential equations, Mittag-Leffler function, Wright 
function. Random time. First passage time. 
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1 Introduction and background 


Fractional Poisson processes and related counting processes are attracting the attention of 
several authors. Most of the recent papers on this topic are centered on certain fractional 
versions (time-fractional, space-fractional, space-time fractional) of the Poisson process, as 
well as some fractional birth processes (see, for instance, the review in [Orsingher (2013')] 
and [Alipour et al (2015)] ). Moreover, [Beghin and Orsingher (2010)| study the properties of 
Poisson-type fractional processes, governed by fractional recursive differential equations, ob¬ 
tained substituting regular derivatives with fractional derivatives. [Mainardi et al (2004)| pro¬ 
vide a generalization of the pure and compound Poisson processes via fractional calculus, by 
resorting to a renewal process-based approach involving waiting time distributions expressed in 
term of the Mittag-Leffler function. A different approach has been developed by [Laskin (2003)] 
and Laskin (2009)| , where a fractional non-Markov Poisson stochastic process based on a 
fractional generalization of the Kolmogorov-Feller equations, and some interesting applica¬ 
tions including a fractional compound Poisson process have been considered. More recently. 
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[Meerschaert et al (2011)] show that a Poisson process, with the time variable replaced by 
an independent inverse stable subordinator, is also a fractional Poisson process. Other re¬ 
cent results on fractional Poisson process can be found in [Gorenflo and Mainardi (2012)| and 
[Gorenflo and Mainardi (2013)| . 


Gounting processes with jumps of amplitude larger than 1 are employed in various applica¬ 
tions, since they are useful to describe simultaneous but independent Poisson streams (see 
[Adelson (1966)| , for instance). The case of fractional compound Poisson processes has been 
investigated by [Scalas (2011)], [Beghin and Macci (2012)] and [Beghin and Macci (2016)] , for 
instance. Moreover, [Beghin and Macci (2014)| consider two fractional versions of nonnegative, 
integer-valued compound Poisson processes, and prove that their probability mass function 
solve certain fractional Kolmogorov forward equations. Gertain fractional growth processes in¬ 
cluding the possibility of jumps of amplitude larger than 1 have been obtained recently through 
the interesting space-fractional Poisson process (cf. [Orsingher and Polito (2012)]) and, more 
generally, through the class of point processes studied in Orsingher and Toaldo (2015)] and 
[Polito and Scalas (2016)'] . The relevance of fractional compound Poisson processes in applica¬ 
tions in ruin theory and their long-range dependence are investigated in [Biard and Saussereau (2014)] 
and [Maheshwari and Vellaisamy (2016)]. 


Following the lines of the papers above, here we analyse a suitable extension of the fractional 
Poisson process, say which performs k kinds of jumps of amplitude 1,2,..., A; with 

rates Ai, A 2 , • • •, Afc respectively. (Throughout the paper we refer to the fractional derivative in 
the Gaputo sense, also known as Dzherbashyan-Caputo fractional derivative). We first obtain 
the moment generating function and the probability law of the process, and discuss its equiva¬ 
lent representation in terms of a subordinator governed by a suitable fractional Gauchy problem. 


Along the same lines as [Beghin and Orsingher (2010)] , in Section [2| we consider the difference- 
differential equations governing the probability mass function of M^(t) and involving the time- 
fractional derivative of order u G (0,1]. The solution of the resulting Gauchy problem repre¬ 
sents the probability distribution of the fractional counting process Hence, we obtain 

E and p^(t) = (t) = k} in terms of a generalized Mittag-Leffler function. We 

also show two useful representations for 

(i) We prove that M^{t) can be expressed as a compound fractional Poisson process. This 
representation is essential to obtain a waiting time distribution. 

(ii) We show that can be regarded as a homogeneous Poisson process with k kinds of 

jumps stopped at a random time. Such random time is the sole component of this subordinat¬ 
ing relationship affected by the fractional derivative, since its distribution is obtained from the 
fundamental solution of a fractional diffusion equation. 


In Section [3] we face the problem of determining certain waiting time and first-passage-time dis¬ 
tributions. Specifically, we evaluate the probability that the first jump of size j, j = 1, 2,... , A:, 
for the process occurs before time t > 0. Interestingly, we prove that the first occurrence 

time of a jump of amplitude j has the same distribution as the waiting time of the first event 
of the classical fractional Poisson process defined with parameter Xj, for j G {1, 2,..., A;}. This 
is an immediate extension of a well-known result. Indeed, for a Poisson process with intensity 
Al -|- A 2 and such that its events are classified as type j via independent Bernoulli trials with 
probability , the first occurrence time of an event of type j is distributed as the interarrival 

time of a Poisson process with intensity Xj, j = 1, 2. In Theorem Id.Il we extend this result to the 
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fractional setting. The remarkable difference is that the exponential density of the interarrival 
times of the Poisson process is replaced by the corresponding density of the fractional Poisson 
process, which depends on the (two-parameter) Mittal-Leffler function. In Section [3] we also 
study, when k = 2, the distribution of the first passage time of M^(t) to a fixed level. We ex¬ 
press it in an integral form which involves the joint distribution of the fractional Poisson process. 

Finally, in Section 0] we obtain a formal expression of the moments of and show that 

both the ratios given by the powers of the fractional Poisson process and of the process M^(t) 
over their means tend to 1 in probability. This result is useful in some applications. In fact, 
from a physical point of view, it means that the distance between the distributions of such 
processes at time t and their equilibrium measures is close to 1 until some deterministic ‘cutoff 
time’ and is close to 0 shortly after. 

In the remaining part of this section we briefly recall some well-known results on the fractional 
Poisson process which will be used throughout the paper. Consider the fractional Poisson 
process 

> 0} , e (0,1], A G (0,oo), (1) 

namely the renewal process with i.i.d. interarrival times distributed according to the fol¬ 
lowing density, for j = 1, 2,... and t G (0, oo) (see [Beghin and Orsingher (2010)| ): 

{t) =¥{% Gdt}/dt = (2) 


where 

CO ^ 

a,l3 ec, Re{a), Re{j3) > 0, x G M 

^ T{ar + I3) 

is the (two-parameter) Mittag-Leffler function. From the Laplace transform 

it follows that the density of the waiting time of the k-th event, possesses the 

Laplace transform 

\k 

^{/-(i);4 = -^-.. 

(s- + A)^ 

Its inverse can be obtained by applying formula (2.5) of [Prabhakar (1971)] , i.e. 

r / \ '1 7 

i?{t-'->E^,,(c.t'>);,} = ^-^.^, (3) 

(where Re{j3) > 0, Re^) > 0, Re{6) > 0 and s > |a;| R<P)). By setting fd = u, '^ = ku, 6 = k 
and cj = — A we have 

/(( (t) = pm G dt} jdt = (4) 

where 

OO / \ ^ 

r!r(«r+ /3) ’ ^ ^ ^^A)Ae-{ld),Re-{l) > 0 (5) 
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is a generalized Mittag-Lefiler function and, as usual, (7)^ = 7(7+1 )... (7 + r— 1), r = 1,2,..., 
(7)0 = 1, is the Pochhammer symbol. 


The corresponding distribution function can be obtained by integrating (a, thus obtaining 
(see Eq. (2.20) of [Beghin and Qrsingher (2010)| ) 

Fl;{t) ^F{n<t}= (6) 

Taking into account ([6]), the probability mass function of the process N^(t) can be easily 
computed as follows (see, also, Eq. (2.21) of [Beghin and Qrsingher (2010)| ): 

P{iV):(t) = n} = F{T^<t< r„+i) = (At")”EC+i+i(-At"). (7) 

Moreover, recalling Eq. (2.29) of [Beghin and Qrsingher (2010)] , we have that the moment 
generating function of the process N^(t), t > 0, can be expressed as 


E 




= E,,i(A(e^-l)t"), 


s E 


( 8 ) 


The mean and the variance of (t) read (see Eqs. (2.7) and (2.8) of [Beghin and Qrsingher (2009)]) 


E[ivj:(t)] = 


xe 


r(i. + i)’ 


Var[iVj:(t)] = 


2 (At" 


(At")^ 


+ 


At" 


T{2u + i) (r(z/ + i))2 r(i/ + i) 


(9) 


In general, the analytical expression for the mth order moment of the fractional Poisson process 
is given by (cf. [Laskin (2009)] , Eq. (40)) 

m 

E[(jvnt)n = E®Em.i)(An'. (lo) 

z=o 

where S^, (m, 1) is the fractional Stirling number defined by Eq. (32) of [Laskin (2009)]. 


2 Fractional counting process 

Let {M^(t);t > 0} be a counting process defined by following rules: 

1. Mi(0) = 0 a.s.; 

2. M^(t) has stationary and independent increments; 

3. F{M^{h) = j} = Xjh + o(/i), for j = 1, 2,..., fc; 

4 . F{M^h) > k} = o{h), 

where /c E N = {1, 2,...} is fixed, and Ai, A2 ,..., Afc >0. From the above assumptions we have 
that the probability distribution pj{t) = E|M^(t) = j}, for j £ Nq = { 0 , 1 , 2 ,...}, satisfies 
the following system of difference-differential equations: 

= '^XrPj-r{t) - {Xi + ... + Xk)pj{t), t>0, (11) 

r=l 

where Pj{t) = 0 for j < 0. 

In this section we examine a fractional extension of {M^{t)',t > 0}. We obtain a proper 
probability distribution and explore the main properties of the corresponding fractional process. 
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2.1 The probability law 

With reference to the fractional derivatives 


dV(t) 

dt^ 


1 f* 
r(i-y) Jo 


(d/ds)f(s) 

{t-sV 


ds, 


0 < V <1, 


f'{t) 2^=1, 


let us now introduce a fractional extension of the process For all fixed u G (0,1] and 

fe G N, let > 0} be a counting process, and assume that the probability distribution 


pJ(f) = P{M"(f)=i}, jGNo 

satisfies the following system of fractional difference-differential equations 

'dpoW 


dt^ 

dp’jjt) 

dt^ 

dp^(f) 


= -Apo(t) 

= Er=l KPj-rit) - ^Pj{t), 
= Er=l ^TP'--rit) - ^pUt), 


for A = Ai -|- A2 + ... + Afc, together with the condition 

"l, j = o 


Pi(0) = 


j = 1,2,... ,fc - 1 

j = A:, A; -Fl,..., 


0, j>i. 


( 12 ) 


(13) 


(14) 


Clearly, when u = 1 the system m identifies with the difference-differential equations of 
process M^{t) given in dill) . Furthermore, when A; = 1 the process identifies with the 

process N^{t) considered in Section [TJ 

Hereafter we will obtain the solution to (I13p - (|14p in terms of the generalized Mittag-LefHer 
function ([5]) and show that it represents a true probability distribution of To this 

purpose we first obtain the moment generating function of in terms of the Mittag-Leffler 

function. 


Proposition 2.1. For all fixed v G (0,1] and A G N, the moment generating function of M‘'(t) 
is given by 

k 

= F;^,i(^Aj (e^'" - 1) t>0, sGM. (15) 

i=i 

Proof. From system (1131) and condition (I14p we have that the probability generating function 
G{z,t) := E satisfies the Cauchy problem 


E 




i=i 

G{z,0) = 1. 

By adopting a Laplace-transform approach we obtain 


^{G{z,t)-s} 


S-+ E •=! A, (l - 


Eq. (jl5l) thus follows recalling Eq. ©• 


□ 
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We remark that the use of the Caputo fractional derivative permits us to avoid fractional initial 
conditions in the previous proof since, in general, 




Z.G (0,1], 


a:=0 

Let us now show that M^(t) can be expressed as a compound fractional Poisson process. 
Proposition 2.2. For all fixed v G (0,1] we have 


M^{t) = 


N’fiit) 

E 

i=l 


Xi 


t > 0, 


(16) 


where N^{t) is a fractional Poisson process, defined as in with intensity A = Ai + A2 + 
.. . + Afc. Moreover, {Xn : n > 1} is a sequence ofi.i.d. random variables, independent of N^{t), 
such that for any n G N 

F{Xn=j} = ^, j = l,2,...,k (17) 

and where both N^(t) and Xn depend on the same parameters Ai, A2, • • •, A^. 

Proof. The moment generating function of Y{t) := Xi, t > 0, can be expressed as 


E 




= E 
= E 


E 


,sY{t) 


(E [e"^']) 


xm 

Nm 


Hence, since 


we have 



i=i 


E 


^sY(t) 


= E 


=iVXWln(iEj=iA, 



Finally, making use of Eq. ([8]) we immediately obtain that the moment generating function of 
Y (t) identifies with the right-hand-side of (fTSll . This completes the proof. □ 


We remark that, due to Proposition 12.21 M^(t) can be regarded as a special case of the pro¬ 
cess defined in Eq. (7) of [Beghin and Macci (2014)] , under a suitable choice of the prob¬ 
ability mass function {qk)k>i the parameter A. Eurthermore, according to Definition 
7.1.1 of [Beghin and Korolev (2002)| , the process M^(t) is a compound Cox process, since 
[Beghin and Orsingher (2010)| show that (t) is a Cox process with a proper directing mea¬ 
sure. Moreover, M'^{t) is a compound fractional process, and thus it is neither Markovian nor 
Levy (cf. jScalas (2011)] ). 

We are now able to obtain the probability mass function (11211 of M'^{t). Indeed, the following 
Proposition holds true. 

Proposition 2.3. The solution pj ft) of the Cauchy problem for j G No, G (0,1] 

and t >0, is given by 


m =E E 

r=0 ai+a 2 +...+ak=r 

ax+2a2+...+kak=j 





r+1 


{-kC). 


(18) 
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(a) (b) 

P)(!) p'jW 



Figure 1: Probability distribution of given in (1181) . for j = 0,1,..., 11, with k = 3, 

V = 0.5, Ai = A2 = A3 = 1, (a) t = 1 and (6) t = 2. The displayed probability mass is (a) 
0.797292 and (6) 0.629278. 


Proof. From (11611 and from a conditioning argument we have 


j 

= P {M" {t)=j} = Y,^{Xi+X2 + ... + Xr= j} P {iVX(t) = r} . 

r=0 

Since Xi,X 2 , ■ ■ ■ ,Xj. are independent and identically distributed (cf. (1171) 1. it follows that 


F{Xi+X2 + ...+Xr 


J} 


E 

OLi+a2+...+ak=r 

ai+2a2+...+fcofe=i 



X 





where the sum is taken in order to consider all the possible ways of having r jumps, with ai 
jumps of size 1, ..., jumps of size k, and such that the total amplitude, i.e. ai+2a2+- • -Pkak, 
equals j. Hence, recalling formula ([7]), the Proposition follows. □ 

Proposition 12.31 is an extension of Proposition 2 of |Di Crescenzo et al (201^ , which is con¬ 
cerning with case k = 2. Some plots of probabilities (I18p are shown in Figure [1] and Figure 

El 

From (jlSp we note that, for v G (0,1], 


p^it) = t>0. 

Moreover, making use of Eqs. ([5]) and (I18p we obtain hereafter the distribution of the process 
in the special case v = 1. 

Corollary 2.4. The probability mass function p]{t), for j G Nq and t >0, is given by 




A“iA“E..A“'= 


r=0 ai+a 2 +...+ai^=r 

ai+2a2+...+kak=j 


oi! 02 ! ... Ofc! 


k 


(19) 
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(a) 


(b) 

p]i’) p]m 



Figure 2 : Probability distribution of given in (fTHIl . for 0 < t < 2 , with A: = 3, Ai = A 2 = 

A 3 = 1, (a) v = 0.5 and ( 6 ) 1 ^ = 1. 


2.2 Equivalent representation 

We will now examine an interesting relationship between the process M^{t) and the process 
In fact, we show that the following representation holds: 

= (t)), 


where (t) is a suitable random process, and thus can be considered as a homogeneous 

Poisson-type counting process with jumps of sizes 1, 2,... , /c stopped at a random time (t)- 
Let us denote by g{z,t) = g2u th® solution of the Cauchy problem 


_ d^g(z,t) 

dT^ dz ^ ’ 


< 9{z,0) = S(z), 

dg(z,t) _ „ 

dt — 


t=o 


t > 0, z G M 
0 < 1/ < 1 


It is well-known that (see [Mainardi (1996a)] and [Mainardi (1996b)] ) 


g2v iz,t) = —W-u,l-v 



t > 0 , 2; G M, 


where 

{x) 

is the Wright function. Let 


00 ^ 

E X 

k\ P {ak + ( 5 ) ’ 


a > —1, /3 > 0, X G K 


g2u {z,t) 


2g2u{z,t), z >0 

0 , 2; < 0 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


be the folded solution to pHll and let (t) be a random process (independent from the process 
(t)) whose transition density F{^i, (t) G dz} jdz is given in (|23p . 

Remark 2.5. It has been proved in [Orsingher and Beghin (2004)] that the solution g2u to (I20p 
can be alternatively expressed as 


g2u (z, t) 


1 


{t — w) {w,\z\) dw, 2; G M, 


2r(l - iz) 


0 























where fu {■,y) is a stable law Siy {y,l3,a) of order with parameters y = 0, (3 = 1 and a = 
(zcosf)i 

Proposition 2.6. The process M'^ (t) and the process (t)) are identically distributed. 

Proof. From (fT^ and (l2Hll we have 


POO 

F {M^ {£^ 2 u (t)) = n} = / pl{z)g 2 u{z,t)dz. 

Jo 

Hence, making use of (I19p and (12211 we get 


’ i^2u (t)) = n} =Y^ 


\oi \a2 \oik 
Ai A 2 •• • A^ 


oi!oo! • • • oiP. 

3=0 01+02+...+Ofe=J 

ai+2a2+...+kak=n 


1 r 
^ Jo 


e z^W-^^i-u (“I:) dz. 


For y = Az, the last expression identifies with (|18p due to the following integral representa¬ 
tion of the generalized Mittag-Leffler function in terms of the Wright function, proposed by 
[Beghin and Orsingher (2010)|: 


Tpk+l 

^v,ku+l 


{-Kn 


1 

k\ 






This completes the proof. 


□ 


Remark 2.7. Since the transition density (I23p coincides with the probability density function 
of the standard inverse z/-stable subordinator E^(t) (see [Meerschaert et al (2011)1 ), the result 
given in Proposition 12.61 can be stated also as follows: The process (t) and the process 
(E^(t)) are identically distributed. 

Remark 2.8. In [Beghin and Orsingher (2010)] Beghin and Orsingher proved an analogous sub¬ 
ordination relationship, i.e. 

JVi(+ 2 +t)), 

where Nf^ft) is the fractional Poisson process defined in ([1]) and f^ 2 u{t) is the random time 
defined above. 

Remark 2.9. By taking u = ^, from Proposition 12.61 we have that (t) and (t)) 

are identically distributed. We note that the random time (t), t > 0, becomes a reflecting 
Brownian motion. Indeed, in this case equation (12011 reduces to the heat equation 


dz ^ 

g{z,0) =d{z), 



t > 0, z G M 


and the solution gi [z, t) is the density of a Brownian motion B (t) ,t > 0, with inhnitesimal 
variance 2. After folding up the solution, we find the following probability mass 


P{M^(^i(t)) 



_2_ 

e 4t 


dz 


= F{M^{\B{t)\) = n}, 


so that (t) is a jump process at a Brownian time. 
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Remark 2.10. It is worth noticing that both the compositions of the fractional Poisson process 
N^{t) defined in ([T]) and of the fractional process M^{t) defined in (|12p with the random time 
^ 2 v {t) yields again fractional processes of different order, i.e. 

iVA(^2. (t)) = Nf{t) and (t)) = 

Taking into account the subordinating relationships examined in Proposition 12.61 and in Re¬ 
mark 12.81 this fact follows immediately from Remark 3.1 of [Kumar et al (2011)] , since, in 
general, the composition of two stable subordinators of indexes /3i and /32 respectively is a 
stable subordinator of index /3i/32- 

Remark 2.11. Bearing in mind Proposition 12.21 setting 

k 

Sr = A- E[xn = r = 1, 2 

j=i 

and recalling ([9]), we can compute more effortlessly the mean and the variance of the process. 
In fact, by Wald’s equation we have 


E [M^(t)] = E[X] • E [iVX (t)] 


Si 

r>TT)’ 


t > 0 . 


Moreover, by the law of total variance we get 


Var[M^(t)] 


Var [X] • E [Nm + (E [X])^ • Var [N}i{t)] 


S2t^ 

r(i. + i) 




t > 0, 


where 

u (r(2z^) ~ i/r2(iy) 

As a consequence it is not hard to show that Var [M’^(t)] — E[M'^(t)] > 0, or, equivalently, 
that the process exhibits overdispersion, since Ziyu) > 0 for all v E (0,1) and Z{1) = 0. 

Finally, we point out that a formal expression for the moments of process is provided in 

Lemma O 



3 Waiting times and first-passage times 


We evaluate the probability distribution function of the waiting time until the first occurrence 
of a jump of size f, f = 1, 2,..., A:, for the process We first observe that the following 

decomposition holds: 

k 

= t>0, 

i=i 


where 


A’XW 

•= j = l,2,...,/c. 


(24) 
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and thus M^{t) counts the number of jumps of amplitude j performed by in (0,t]. 

Furthermore, we introduce the random variables 

Hj := inf {s > 0 : MJ(s) = l} and ~ Geo > i = Ij 2,..., fc. 

In other words, Hj represents the first occurrence time of a jump of amplitude j for process 
whereas is a geometric random variable with parameter ^ that describes the order 
of the first jump of amplitude j in the sequence of jumps of We prove that Hj is 

distributed as the waiting time of the first event of the fractional Poisson process defined in © 
with parameter Xj. Indeed, the following result holds. 

Theorem 3.1. Let j E {1, 2,..., A;}. Then 


^{Hj <t} = XjHE { — XjH) , t > 0. 

Proof. By conditioning on G^, for t > 0, due to Eqs. (US]) and dSj) we have 
F{Hj <t} = Egj [P {Hj < t I G^}] 

H-cxd 

= <t \ G^ = n} P [G^ = n} 

n—1 


(25) 


n=l 

H-oo 




n=l 

H-cxd 


21 

A 




n=l 


A 


n—1 


= E A-c- 1 - ^ 

n=0 ^ ' 

By using formula (2.3.1) of [Mathai and Haubold (2008)] , i.e. 

lo ~ ^ (^) ’ 

(where Re{a) > 0,Re{/3) > 0 and Re{'y) > 0) for a = (3 = n, ■j = nn + 1, 6 = n + 1 and 
z = —AG, we get 


^{Hj<t} = 


XiG 

W) 


l-J / a"-(l-ur‘(-AtV),!!, 

n=0 ^ /Jo 


Xjt 

r>) 


-j-CXD 


u rl _ 

/ {i-’-r'E 

Jo 


n=0 


AG{l---^)u^ E::+l^,{-AGundu. 


Due to formula (2.30) of [Beghin and Orsingher (2010)] , i.e. 

-j-CXD 

{XwGTe:^^^, i-XG) = (A {w - 1) G) (hi < 1, t > 0), 


n=0 


we have 


F{Hj <t} = 


X,G 


(1 — n)^ ^ E,y^i{—XjGu'^)du 


T{v) 
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By making use of formula (2.2.14) of [Mathai and Haubold (2008)| , i.e. 

[ (1 - zy~^ Ea,i 3 (xz^^) dz = T {a) (x ), 

Jo 

(where a > 0; /3,a G C] Re{/3) > 0 and Re{a) > 0), for a = a = u, /3 = 1 and x = —XjE, we 
get 

P {Hj <t} = XjEE^^y+i (-XjE), t>0. 

Therefore Hj is distributed as the waiting time of the first event of the fractional Poisson 
process dehned in ([T|) (cf. Q). □ 

The result shown in Theorem ED is an immediate extension of the well-known result for the 
Poisson process, i.e. for i/ = 1, by which Hj is exponentially distributed with parameter Xj. 

We will now be concerned with the distribution of the hrst passage time to a fixed level for the 
process denoted as 


Tn = inf {s > 0 : {s) = n} , n G N. (26) 

The following result is concerning the case k = 2, i.e. when the process performs jumps 

of sizes 1 and 2. 


Theorem 3.2. The cumulative distribution function of the first passage time when k = 2 
reads 


+ 00 h 


i*K<i)=E E E 




z 

n — i 


3 - * 


h - n - j + ij V ■^1 + ^2 


Ai 


2j-h 


X 2 


Ai -|- A2 


h-j 


f p (t) = j, (s) = z} ds, t > 0. 

Jo 


(27) 


Proof. Since the process M'^{t) performs jumps of size 1 and 2, and has non-independent 
increments, the computation of the cumulative distribution function of the hrst passage time 
(1261) can be carried out as follows: 


- 1-00 

P{rn<t}=V/ r (t) = h, (s) = n} ds 

h=J^ 

- 1-00 h j 

= E E E / P (t) = h, (s) = n I (t) = J, (s) = i} 
X p{ivj:^+,, it)= i,ivj:i+A, (^) = *} ds. 
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Making use of Proposition 12.21 we have: 


’ {Tn <t} = 


+00 h j ( j i 'j 

E E E/i* E^'' = '>'E^' = " 

/i=nj = |-il i=l lr=l 1=1 ) 


X P {t) = 3 , (s) = i] ds 

+00 h j ( i j 'j 

p|^Xz=n, ^ Xr = h-n\ 
h=nj=|'.|'| j=l I/=! r=i+l ) 

X P it) = j, is) = i} ds 

+ 0 O h j ( i ^ ^ 

^EEE/r E^. = 4 He 

h=n i=l I/=! ) l^r=«+l 

X P {t) = J, {s) = i] ds 


Xr = h — n 


y y y ^ 

u u \n - ij \Xi + ^2 

h=nj=^h^ *=1 


A 2 


J - i 


h - n - j + ij V -^1 + -^2 


Ai + A 2 

A 2 


^ 2j—2i-\-n—h / \ \ h—n—j+i 


Ai + A 2 


X / ^{Nl^+xAt)=j^Nl^^^^{s)=i]ds, 
Jo 


this giving Eq. (127)1 . 


□ 


To the best of our knowledge, the bivariate distribution shown in the right-hand-side of (I27h . 
i-*' f {JVJ.+A. (t) = j}, cannot be expressed in a closed form. [Orsingher and Polito (2013)| 

derived an expression in terms of Prabhakar integrals, i.e.: 


PjiVAVA, (4=bA^A4+A, {t)=j] 

~ T A 2 ) 

xy"-iE,,,(-(Ai + A2)2/4)(^))(4, 

where 

(a^) = J (^ - “ 4^) </> (4 dt 

is the Prabhakar integral (see [Prabhakar (1971)] for details). Politi et al. jPoliti et al (2011)] , 
instead, evaluate the joint probability given in ()27p by introducing the random variable Yi 
which denotes the residual lifetime at s (that is the time to the next epoch) conditional on 

('®) ~ b i-®- [t ~ s I (s) = i] whose cumulative distribution function is 

denoted by Fy^iu)- Therefore, 

P {a^a4+a, (^) = b it) = j} =P (t) - (4 = i - * I a'a.+a, (^) = 4 

xP{1V,4+A2 (4 = ^}, 
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where 


P {iVj:,+A, it) - +A, is)=j-i\ (s) = i} 

^ f Io~" {^Ai+A 2 it-s-y)=j -i-l} dFYi (y), if j - i > 1, 
|l-Fy,(i-s), ifj-i = 0. 

It is meaningful to stress that when k = 2 the passage of M^{t) to a level n is not sure. In fact, 
the process can cross state n without visiting it due to the effect of a jump having size 2. 

4 Convergence results 

For the processes N^(t) and defined respectively in ([1]) and in (I12p . we now focus on a 

property related to their asymptotic behavior as the relevant parameters grow larger. 

Proposition 4.1. Let n G (0,1]. Then for a fixed t > 0 we have 

A^A (^) Prob ^ ^ 

E[N^{t)] A^+oo 

Proof. We study the convergence in mean of the random variable to 1- Due to the 

triangular inequality we have 


E 


E [iV-(t)] 


< 2. 


Therefore, we can apply the dominated convergence theorem and calculate the following limit: 


lim E 

A—^~l“Oo 


A^A"(t) 


E [ivj:(t)] 


-1 


H-OO 

: lim 
A—>+oo 

j=0 




- 1 


r(i.+i) 




(28) 


Taking account of the behavior of the generalized Mittag-Leffler function for large z (see 
[Saxena et al (2004)] for details), i.e.: 


~ o (|; 


I -(5 


\z\ > 1, 


we can conclude that limit (I28p equals 0. This fact proves the result since convergence in mean 
implies convergence in probability. □ 


The previous result can be extended to a more general setting. Recalling the expression (llUp 
for the moments of (t), the proof of the next proposition is similar to that of Proposition 
O and thus is omitted. 

Proposition 4.2. Let u G (0,1] and r G N. Then, for a fixed t > 0, 


[A^a"( 0]' 


Prob 


E{[jvut)]q 


A—^-|-oo 


> 1. 


In order to prove an analogous result for M^{t), in the following lemma we give a formal 
expression for the moments of such a process. 
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Lemma 4.1. The order moment of the process t > 0, reads 

tru 


E{iM''{t)r} 


E 


r (ri' + 1) 
r=0 ii+...+i]^=r 


X E 

ni+...+nfc=m 


m 

Ti \, • • • , Tl}^ 


“ firik 

dgnr\ ) ds’^k'- > 


(29) 


s=0 


Proof. By applying Hoppe’s formula in order to evaluate the derivatives of the moment gen¬ 
erating function of the process cf. (fT^ . we have 


^ l2=V'' X iei^-Dt^ 

E{[M’^(t)]”"} = ^- E,.i 

r =0 


ri 


Am,ri 

Vi=i 


s=0 


where 


^m,r I 


EAi(e^”-l)i‘' =E 


-1 = 1 


h=0 


ds'^ 


( k \ 1 —h 

- E A; (<"■ - 1) j 

k \ ^ 

E Aj («'■■ - 1) C . 

/ 


Finally, after using rather cumbersome algebra, we obtain (1291) . 


□ 


It is now immediate to verify the following result for M^(t). 


Proposition 4.3. Let ix G (0,1] and m G N. Then, for i G {1, 2,... ,k} and for a fixed t > 0, we 
have 


[M-{t)r Prob . 


Proof. By virtue of (1291) . convergence in probability can be obtained by proving convergence 
in mean, as in Proposition 14.11 □ 


The results presented in this section deserve interest in some physical contexts. We recall that 
a family of random variables exhibits cut-off behavior at mean times if (see, for instance, 
Definition 1 of jBarrera (2009)|) 

Prob ^ ^ 

E [17(^)] A^+oo^ 

Hence, Propositions 14.1114.21 and 14.31 show that the processes [N^ (i)]™ and [M^ m G N, 

exhibit cut-off behavior at mean times with respect to the relevant parameters or, roughly 
speaking, that they somehow converge very abruptly to equilibrium. 


We finally remark that in this context the sufficient condition given in Proposition 1 of [Barrera (2009)1 
is not useful to prove Proposition 14.11 since such condition holds only when n = 1. 
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